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THE  STEADY  MOTION  OF  A SYMMETRICAL 
OBSTACLE  ALONG  THE  AXIS  OF  A 
ROTATING  FLUID 

ABSTRACT 

An  investigation  by  G.  I.  Taylor  of  the  steady  motion  of  an  obstacle  along 
the  axis  of  a rotating  fluid,  is  extended  in  the  following  paper.  It  is  3hown 
that  Taylor's  particular  solution  is  just  one  of  an  infinity  of  functions 
comprising  the  general  solution. 

The  theory  is  applied  to  motions  in  a rotating  cylinder  of  fluid.  A 
critical  Rossby  number  is  derived,  below  which  the  flow  around  the  obstacle  is 
wave-like.  When  the  Rossby  number  is  greater  than  the  critical  value,  the  flow 
consists  only  of  a local  perturbation  that  dies  out  rapidly  on  both  sides  of  the 
obstacle.  Various  other  critical  numbers  exist,  below  which  additional  mode 3 
of  oscillation  become  dynamically  possible. 

An  experiment  was  designed  to  test  the  theoretical  results  of  this  paper. 

An  obstacle  was  moved  along  the  axis  of  a long  cylinder  of  rotating  water.  The 
resulting  flow  patterns  were  observed  visually  and  photographically.  The  three- 
dimensional  wave  motions  which  occurred  in  the  experiment  were  unquestionably 
the  3ame  as  those  in  the  theoretical  solution. 


?:•£  STEADY  ''(TIC!!  OF  A SY! ’"ETHICAL  C35TACI.E  ALONG 
THE  AXIS  C?  A HCTATIJC  FLUID 


1.  Introduction 

In  a aeries  of  naper3  written  a number  of  years  ago,  Taylor  (1Q20,  1922, 
1923),  Proudnan  (1916),  and  Grace  (1922,  1923,  1924,  1926)  investigated  the 
problem  of  the  motion  of  a perfect,  rotating,  fluid  around  obstacles  (especially 
soheres).  This  work  has  acquired  renewed  interest  in  the  recent  work  of  "organ 
(1951)  end  Stewartson  (1952). 

With  the  exception  of  Taylor,  the  problem  is  of  academic  interest  alone 
to  the  above  authors.  Indeed , it  renreser.t3  a fascinating  challenge  for  the 
extension  of  the  investigations  of  classical  hydrodynamics  to  an  area  that  is 
almost  unexplored.  If  the  oroblem  had  no  other  importance,  the  present  paper 
should,  no  doubt,  have  been  submitted  to  the  Proceedings  of  the  Royal  Society 
of  London  in  which  almost  all  of  the  above  researches  have  aDDeared.  It  is 
almost  obvious  that  this  is  not,  or  should  not  be,  the  ca3e.  Certainly  all 
meteorologists  are  agreed  that  the  moderate  and  large-scale  motions  of  the 
atmosphere  are  profoundly  influenced  by  the  vorticity  of  the  earth's  rotation. 

It  would  seem  to  be  evident,  then,  that  efforts  to  understand  these  phenomena 
should  begin  with  studies  of  flow  of  a rotating  fluid  around  elementary  bodies. 
The  motion  of  an  irrotational  liquid  oand  a sphere  is  discussed  in  an  early 
chapter  of  Lamb's  Hydrodynamics,  and  Stokes'  solution  for  viscous  flow  around  a 
sphere  vas  one  of  the  first  investigations  based  on  the  Kavier-Stokes  equations 
of  motion. 

It  is  unfortunate  that  basic  investigations  of  a rotating  fluid  are  infre- 
quent in  meteorological  literature.  The  resulting  lack  of  knowledge  of  the  moat 
fundamental  effects  of  rotation  on  fluid  motions  (.aa  been  a distinct  handicap  in 
current  efforts  to  produce  laboratory  models  of  atmospheric  phenomena.  These 
difficulties  are  discussed  in  a series  of  papers  (Long,  1951,  Long,  1952a,  Fultz 
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One  of  the  most  interesting  of  these  basic  problems  to  a me teorologir-t  is  the 
motion  of  an  obstacle  in  a rotating  fluid  ir.  a direction  perpendicular  to  the  axis 
of  rotation.  This  motion  would  nave  analogies  to  the  effect  of  mountain  barriers 
on  the  ever-present  zonal  currents  In  the  earth's  atmosphere. 

Efforts  in  this  direction  have  not  proved  to  be  particularly  successful, 
however,  since  the  analysis  is  very  complicated.  The  steady  motion  of  a sohere 
along  the  axis  of  a rotating  fluic,  as  analysed  by  Taylor  (1922),  is  simpler  and 
his  re3ult3  were  the  only  ones  that  d’.d  riot  require  the  assumption  of  infinitesimally 
small  motions.  A3  shown  below,  Taylor's  solution  is  only  or.e  of  an  infinite  sum  of 
functions  which  comprise  the  general  solution,  and  he  was  unable,  therefore,  to 
obtain  a physically  realizable  flow  pattern.  In  the  present  paper,  section  2 is 
deveted  tc  a further  theoretical  discussion  of  this  problem.  An  experiment  to  test 
the  theory  is  described  in  section  3. 

2.  General  solution 

Figure  1 is  a sketch  of  the  coordinate  system  used  in  the  theoretical  develop- 
ment of  the  flow  around  a symmetrical  obstacle  moving  steadily  along  the  axi3  of  a 
rotating  fluid.  The  x-axis  is  along  the  axis  of  rotation,  and  velocities  in  this 
direction  are  denoted  by  u.  The  component  of  the  velocity  outward  i3  v,  and  w is 
the  component  of  the  absolute  velocity  tangent  to  circles  with  centers  on  the 
x-axia.  The  center  of  the  coordinate  system  is  assumed  to  move  with  the  obstacle 
at  a aneed  uQ,  so  that  the  flew  with  respect  to  this  frame  i3  3teady.  Obviously, 
thi3  requires  the  additional  assumption  that  the  fluid  extends  to  infinity  in  both 
directions  along  the  &xi3  of  rotation. 

The  velocity  components  are 

dx  d p dA 

u,dt  ’ v = dt  ’ t ’ 


and  the  individual  tine  derivative  is 


d d d 

— * U + V 

dt  dp 


( 2 ) 


In  (2)  the  nartial  derivative  with  resDect  tc  tine  is  missing  because  of  the  issued 
steady  state.  The  derivative  with  resoect  to  ^ is  also  zero  since  the  motion  will 
be  symmetrical  about  the  x-exis^.  The  equations  of  motion  and  continuity  for  an 
incompressible,  homogeneous  fluid  are 


dv  w*  d , p 
= -—  (—  +X), 

(3) 

dt  p dp  q 

du  dp 

;r-77<rx>’ 

(4) 

dw  w v 

— ♦ — = 0 
dt  p 

(5) 

du  dv  V 
T"h  — + — * 0 , 
tf*  dp  p 

(6) 

where  q is  the  (constant)  density,  and  x *s  the  potential  of  the  external  forces 
(gravity).  The  s-quation  of  continuity,  (6),  reveals  that  we  may  express  the  velocity 
comDonents,  u,  v,  by  a stream  function,  *Jv(x,  P ) » as  follows: 

U * ‘ , V* 1 v 7) 

p dp  p dx 

Equations  (2)  and  (7)  show  that,  if  a quantity  i3  individually  conserved  during  the 
motion,  it  is  then  e function  only  of  the  stream  function,  ^ • Combining  terms 


1 This  may  tw  inferred  from  the  fact  that  the  equations  cf  motion  and  continuity, 
and  the  kinematic  boundary  conditions  are  invariant  for  transformations  of  ^ 
to  ^ » const.  This  apDlication  of  group  theory  (Birkhoff,  1950)  is  almost 

trivial  in  this  case. 
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in  (5)  we  find  that 


wp=gty) 


(?) 


We  determine  the  unknown  function,  g(  \fr  ),  by  assuming  that  at  a sufficiently 
great  distance  ahead  of  the  obstacle  the  perturbed  motion  vanishes.  This  requires 

p 

that  g(  \jf  ) » 0.  p_  , where  p is  the  undisturbed  distance  of  the  stream  line 

2 

from  the  axis  of  rotation.  Similarly,  from  (7),  ■=  -u  p /Z  in  a Lagrangian 

o o 

sense.  Combining  these  results, 

<r*  Z£l 


w =' 


(9) 


If  we  define  the  following  quantity: 

e dv  d u 
” dx  dp’ 


(10) 


and  eliminate  the  right  hand  sides  of  (3)  and  (4)  by  cross  differentiation,  we 
obtain 


dt  p P 


(11) 


or 


4*4**W!. 


(12) 


Evaluating  H(  ^ ) in  p.  manner  simile.'  to  the  determination  of  g(  ^ ) in  (8), 

-+o-2(^  + ^2)  = ° * (13) 

p r 2 

The  vorticity  component,  f , may  be  expressed  in  terms  of  <j/  by  using  (7)  and  (10). 
The  following  two  equations  then  replace  the  original  set  (3) -(6): 

p»‘*« r'k  (H) 


dx*  dp*  p dp  * 2 


(15) 


The  general  solution  is 


2 

~p~~  = ' ^ +^|{Akcosks  + Bksinkx}J,[(cr2-|<2j2/)] 


(16) 


In  this  equation  J,  is  a Bessel  function  of  the  first  kind.  The  Neumann  function 
becomes  infinite  on  the  axis  of  rotation  ( p = 0)  and  hence  is  omitted. 

Equation  (15)  may  also  be  expressed  in  snnerical  coordinates  by  the  relations 
(see  fig  1) , 

x=  rcos 9,  p=rs\nB  . (17) 

The  corresponding  form  of  the  solution  is 

2 j » gQ 

'i'  =-Us’r  *,n + (xr)*T  sine  P^(cos0) [a„ JnJcrr)  -*■  Bn  Jj^^(crr)j  . (1^) 


The  Bessel  functions  are  those  of  the  half-odd  order  and,  as  is  well  known,  are 

expressible  in  terms  of  elementary  functions  (Margenau  and  Murphy,  1943). 

1 

(cos  9 ) is  an  associated  Legendre  polynomial.  The  particular  solution  for 
n = 1 is  the  solution  obtained  by  Taylor, 

Taylor  experienced  a difficulty  in  applying  his  one  solution  to  the  case  of 
the  flow  around  a sphere  in  a rotating  fluid  extending  to  infinity  in  all  directions. 
The  kinematic  condition  at  the  obstacle  and  the  condition  that  the  perturbation 
vanish  at  infinity  determined  only  one  of  the  two  constants  at  his  disposal,  and 
he  was  left  with  an  indeterminate  problem.  The  fact  that  the  present  paper  adds 
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a new  infinity  of  solutions  would  seem  to  be  a further  complication.  Experience 
with  other  problems  of  a similar  nature  indicates  that  this  is  not  the  case 
(Rayleigh  1383,  Lyra  1943,  Long  1952b).  In  these  references  fluid  systems  a1"" 
considered  in  which  obstacles  set  up  waves  in  a current  passing  over  or  around 
them.  In  all  cases  the  same  indeterminacy  is  found.  The  first  investigation  of 
this  phenomenon  was  made  by  Rayleigh,  who  discovered  that  the  introduction  of 
viscosity  in  the  equations  of  motion  removed  the  indeterminacy.  Then,  letting 
the  frictional  coefficient  tend  to  zero,  the  solution  is  unique  and  involves  a 
complete  damping  of  the  upstream  waves.  The  resulting  asymptotic  solution, 
called  the  "practical  solution"  by  Lamb  (1932) , is  the  one  that  is  realized  in 

actual  exoeriments  (Bakhmeteff,  1932).  Applying  this  to  (18),  all  the  remaining 

v 

constants  would  be  required  to  annul  the  "upstream"  oscillations.  The  experimental 
work  described  in  this  paper  verifies  that  no  waves  occur  in  advance  of  the  obstacle 
(see  fig  9). 

Although  (16)  or  (18)  represents  the  general  solution  for  the  motion  considered 
in  this  naper,  there  are  great  difficulties  in  finding  a solution  that  satisfies 
all  the  conditions  in  a given  special  case.  Even  for  a sphere  moving  in  an  unlimited 
fluid  the  problem  of  choosing  the  coefficients  of  the  Bessel  functions  in  (18)  so 
a3  to  annul  the ■ "upstream"  waves,  seems  prohibitively  complicated.  If  the  cylinder 
has  a finite  radius,  the  difficulties  are  even  greater.  Yet  the  complete  solution 
for  a simple  type  of  obstacle  is  greatly  to  be  desired  in  order  to  ascertain  its 
behavior  as  the  parameter  er  tends  to  infinity  (i.e.  uq  tends  to  zero).  Taylor's 
experiments  and  those  reported  in  this  pater  show  that  an  obstacle,  moving  very 
slowly  along  the  axle,  pushes  a column  of  fluid  ahead  of  it  with  little  disturbance 
in  the  fluid  outside.  This  is  one  of  the  neculiarities  in  the  behavior  of  a rotating 
liquid  and  should  be  carefully  studied. 

A fairly  good  idea  of  the  motion  in  a cylinder  of  radius  b may  be  inferred  from 
(16).  The  kinematic  condition  at  the  outer  wall  requires  that  the  Bessel  function 
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In  terms  of  this  parameter  the  first  mode  of  oscillation  has  a length 


X = 


2irb 

(R*-|4.68)* 


(23) 


It  is  not  surprising  that  for  large  Rossby  numbers  the  obstacle  should  nroduce 
only  a local  perturbation.  An  infinitely  large  SQ  corresponds  to  a vanishing,  basic 
angular  velocity.  The  motion,  therefore,  should  accroach  potential  flow  as  R^— 00  . 
It  is  obvious  from  the  differential  equation  (15)  that  it  does  so. 

The  waves  that  occur  when  Rq  < .261  are  the  subject  of  the  experimental  study 
described  in  the  next  section,  '/ithout  solving  for  the  motion  around  a given 
obstacle  we  cannot  tell  how  the  wave s combine  to  give  the  total  wave  pattern.  If" 
n^  = 1,  however,  only  one  mode  can  exist,  and,  sufficiently  far  behind  the  obstacle, 
this  wave  should  re  ore  sent  the  only  sensible  disturbed  motion.  -Using  the  second 
zero  of  the  function  J-^,  the  one-wave  system  should  occur  in  the  range, 

.142  < R < .261.  (24) 

o 

A plot  of  the  wave  length  vs.  is  shown  in  fig  2. 

The  above  discussion  has  referred  only  to  the  component  of  the  velocity  in 
planes  through  the  axis  of  the  cylinder.  The  component,  w,  normal  to  these  planes 
is  given  analytically  in  (14)  or  (16).  Physically  the  meaning  of  (14)  is  that 
each  fluid  ring  (centered  on  the  x-axis)  conserves  its  angular  momentum  during  the 
motion.  Therefore,  the  ring  slov3  down  when  expanding,  and  speeds  up  when  con- 
tracting, i.e.  when  v is  positive  or  negative  respectively.  The  result  is  a 
three-dimensional  oscillation  which  has  a phase  difference  of  one-quarter  wave 
length  when  projected  on  Dianes  at  right  angles. 

3.  Experiments 

Taylor's  Investigation  of  the  problem  of  this  paper  included  an  experiment 
vrith  a light  sphere  pulled  along  the  axis  of  a rotating  cylinder  of  water.  His 
results  were  very  meager,  however.  No  definite  wave  motion  was  seen,  probably 
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because  re  trued  such  a 3mell  obstacle  th.it  the  vaves  vere  not  excited  by  the  forced 
rer  tur'rat'  or.,  in  view  of  this  deficiency,  and  of  the  new  theoretical  findings  in 
section  2,  it  r.~<-nied  ndv>  sable  tc  n&kt  a further  experimental  study. 

The  equipment  construcT.eJ  for  iriis  work  is  shown  in  fie  3.  The  long  cylinder 
is  90  m in  length  and  has  a:i  inr.»r  :ianeter  of  29  cs.  Trie  obstacle,  shown  in 
position  at  the  axis,  is  a sphere  with  an  attached  cone  to  form  a tapered  trailing 
edge.  The  pur  nose  of  this  is  to  minimize  the  effect  of  boundary  layer  separation. 
The  spherical  head  has  a radius  of  3.8  cm  and  the  entire  obstacle  has  nr.  overall 
length  of  11.2  cm.  The  cylinder  is  mounted  on  a turnl-able  and  rotated  at  a 
uniform  speed  (usually  10  rpm) . "he  water  assumes  the  angular  velocity  of  the 
vessel  within  a reriod  of  15-20  minutes.  The  obstacle  is  lowered  by  a motor  drive 
into  the  water  and  along  the  axis  of  the  cylinder.  A scale,  parallel  to  this  axis, 
is  placed  just  outside  of  the  cylinder  to  facilitate  the  measurement  of  wave 
lengths.  In  view  of  the  dimensions  used  in  fig  2,  this  length  scale  is  in  units 
of  b,  the  inside  radius  of  the  vessel  (14.5  cm). 

The  problem  of  observing  the  flow  in  this  experiment  i3  complicated  by  the 
basic  rotation  and  the  fact  tiiat  the  motions  are  three-dimensional . The  most 
satisfactory  tracer  vas  s ourple  dye  trail  left  by  a dissolving  pellet  of 
Potassium  permanganate  dropned  into  the  water  Just  before  the  descent  of  the 
obstacle  along  the  axis.  These  pellets,  when  drooped  in  a cluster,  leave  a 
group  of  dye  filaments  which  remain  parallel  to  the  axis  of  rotation  until 
deformed  by  the  passage  of  the  obstacle.  Since  the  filaments  are  material  lines 
they  are  the  stream  lines  of  the  motion  relative  to  a coordinate  system  moving 
with,  the  obstacle  and  rotating  with  the  cylinder,  rhotograohing  these  lines 
necessarily  projects  the  Pattern  on  the  Plane  of  the  camera  but  the  wave  length 
is  preserved  by  this  Projection.  The  three-dimensionality  of  the  motion  is 
recorded  by  taking  simultaneous  pictures  with  two  cameras  at  9C  degrees.  It  did 
not  seem  to  be  necessary,  however,  to  include  both  photographs  in  the  figures  of 
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I — Experimental  equipment.  The  long,  glace 
cylinder  is  filled  with  water  and  rotated 
at  a constant  angular  velocity.  The  obstacle 
i3  lowered  at  a uniforn  speed  along  the  exis 
of  rotation. 


The  behavior  of  the  dye  oelT ets , when  dropping  through  the  water,  is  somewhat 
erratic  and  there  is  considerable  difficulty  in  following  a given  filament  in  a 
photograph  over  an  entire  wave  length.  This  lends  some  uncertainty  to  the 
determination  of  the  wave  length  in  most  of  the  photographs. 

Visual  observations  of  the  motion  around  the  obstacle  confirmed  immediately 
the  preliminary  theoretical  deductions  in  section  2 of  this  paper.  As  the  obstacle 
approached  a fluid  ring,  centered  on  the  axis  of  rotation,  it  expanded  to  permit 
its  passage  around  the  barrier.  The  requirement  for  conservation  of  angular 
momentum  caused  it  to  lose  some  of  its  basic  angular  velocity  as  shown  by  the 
twisting  of  the  dye  filaments  in  a direction  opposite  to  the  basic  rotation.  In 
the  lee  of  the  obstacle  the  fluid  ring  performed  a series  of  oscillations,  expand- 
ing and  contracting  and  losing  and  gaining  some  of  its  spin  around  its  axis. 

Two  complications  have  prevented  a detailed  comparison  of  observed  and 
theoretical  wave  lengths.  In  the  first  place,  the  theoretical  findings  were 
based  on  the  assumption  of  an  infinitely  long  cylinder.  In  the  experiment  the 
wave  length  is  an  appreciable  fraction  of  the  cylinder  length.  It  appears  from 
the  observations  that  this  causes  a marked  shortening  of  the  wave  nearest  the  free 
surface.  Secondly,  in  the  vicinity  of  the  obstacle  a local,  perturbation  is  added 
to  the  free  wave  motion.  This  causes  the  first  measured  wave  length  in  the  lee 
of  the  barrier  to  be  somewhat  longer  than  the  theoretical  length.  These  two 
factors  make  the  measurement  of  the  free  wave  length  a subjective  matter, 
especially  for  the  longer  waves.  The  experiments  leave  little  doubt,  however, 
that  the  observed  waves  are  the  same  as  the  theoretical  oscillations.  This  is 
shown  rather  convincingly  in  the  photographs. 

figure  4 is  a photograph  of  the  wave  motion  at  a low  Rossby  number  (0.083). 

In  this  experiment,  as  well  as  in  figs  5 and  6,  RQ  is  in  a range  where  more  than 
one  mode  of  oscillation  is  possible.  Nevertheless,  the  wave  lengths  approach 
those  given  by  equation  (23).  If  the  other  nodes  exist  they  must  have  very  small 
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amplitudes.  Figures  7 and  8 have  higher  Ro3sby  numbers  and  the  waves  are  longer. 
Wave  motion  is  still  evident  at  Rossby  numbers  above  0.171  but  they  are  so  long 
that  it  is  impossible  to  find  one  that  is  not  affected  by  the  local  perturbation 
at  the  obstacle  or  by  the  presence  of  the  free  surface. 

Figure  9 is  a photograph  showing  the  motion  well  ahead  of  the  obstacle. 

As  in  all  experiments,  there  is  no  wave  motion  "upstream." 

Finally  brief  mention  may  be  made  of  the  flow  at  very  low  Rossby  numbers. 

If  Rq  is  less  than  about  0.03  the  experiments  indicate  that  an  obstacle  pushes 
a column  cf  fluid  ahead  of  it  as  it  moves  along  the  axis.  The  radius  of  this 
column  is  indefinite  but  is  somewhat  smaller  than  the  radius  of  the  obstacle  so 
that  some  fluid  still  flows  around  the  barrier.  When  the  -fluid  rings  converge  on 
the  lee  side  a rather  strong  cyclonic  vortex  is  formed  there.  Waves  of  small 
amplitude  and  very  short  wave  length  can  still  be  seen  on  the  "downstream"  side 


of  the  obstacle 
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Fig.  4 — Motion  around  an  obstacle  moving  along  the 
axis  of  rotation  of  a long,  rotating  cylinder 
of  water.  The  Hossby  numoer  of  vhio  experiment 
i 3 0.r''83.  If  the  first  mode  of  oscillation  is 
the  only  one  present,  the  theoretical  length  of 
the  vave3  is  0.31b.  The  waves  between  the  srrows 
have  successive  iength3  of,  aDoroximately,  C.69b, 
0.63b,  0.6lb,  0.57b,  0.53b. 
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Fig.  5 — Wave  motion  around  an  obstacle  ip  » 1 ong, 

rotating  cylinder.  The  Rossby  number  in  thi3 
experiment  is  0.101.  The  theoretical  wave 
length  of  the  first  mode  of  oscillation  is 
0.69b.  The  waves  between  the  arrows  have 
successive  lengths  of  0.70b,  0.69b,  0.b9b. 
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Fig.  6— ’'lave  motion  around  cn  obstacle  in  a long, 

rotating  cylinder.  The.  Rossby  number  in  this 
experiment  i3  0.138.  ne  theoretical  wave 
length  of  the  first  node  of  oscillation  is 
1.02b.  The  waves  between  the  arrows  have 
successive  lengths  o*  ..25b  and  1 .02b. 


Fig.  7 — Wave  motion  around  an  obstacle  in  a long, 

rotating ^cylinder.  The  Rossby  number  in  this 
experiment  is  0.153.  The  theoretical  wave 
length  is  1.19b.  The  waves  between  t}ie  arrows 
have  successive  lengths  of  1.34b  and  1.17b. 
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Fig.  8 — Wave  motion  in  a long,  rotating  cylinder 
of  water.  The  Rossby  number  in  this  exper- 
iment is  0.171.  'he  theoretical  wave  length 
is  1.42h.  The  length  of  the  wave  between  the 
arrows  is  clo.e  to  this  value.  The  shorten- 
ing of  the  waves  near  the  free  surface  is 
quite  evi . ent. 
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Fig.  9 — Motion  around  an  obstacle  in  a long,  rotating 
cylinder  of  water.  This  photograph  3hov3  that 
no  wave  motion  occurs  ahead  of  the  obstacle. 
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